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Abstract 



We present results of a large-scale simulation for the flavor non-singlet light 
hadron spectrum in quenched lattice QCD with the Wilson quark action. 
Hadron masses are calculated at four values of lattice spacing in the range 
a ~ 0.1 - 0.05 fm on lattices with a physical extent of 3 fm at five quark 
masses corresponding to niT^/mp 0.75 - 0.4. The calculated spectrum in the 
continuum limit shows a systematic deviation from experiment, though the 
magnitude of deviation is contained within 11%. Results for decay constants 
and light quark masses are also reported. 
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Typeset using REVT^ 



Deriving the light hadron spectrum is a key step of lattice QCD in order to gain a 
convincing demonstration on the validity of QCD for strong interactions in the low-energy 
non-perturbative domain. Under the present limitation of computer power, the first possible 
step toward this goal is to establish the spectrum in quenched QCD, in which the reaction 
of dynamical sea quarks is switched off. The most extensive effort in this step, amongst 
numerous studies to date has been made by Butler et al. using the Wilson quark 
action. Their results, extrapolated to the continuum limit and corrected for finite lattice 
size effects, show an agreement with the experimental spectrum within 6% for 7 hadron 
masses calculated as compared to the estimated error of 1 to 9% depending on particles. 
A similar attempt with the Kogut-Susskind quark action ||^ has reported a nucleon mass 
higher than experiment by 3% with an error of 4%. 

In this article we report on our attempt at a quenched calculation with a precision 
significantly improved over those of previous studies, thereby establishing the quenched 
spectrum and simultaneously exploring the limitation of the quenched approximation. This 
effort has been made possible by the CP-PACS computer, a massively parallel system with a 
peak speed of 614Gflops developed at the University of Tsukuba ||^. For preliminary reports 
of this work, see Ref. 0. 

We carry out our study with the plaquette action for gluons and the Wilson action for 
quarks. Parameters employed in our simulations are summarized in Table |. Four values 
of the coupling constant (3 = Q/g"^ covering the range of lattice spacing a ^ 0.1 - 0.05 fm 
are chosen, closer toward the continuum limit than a ~ 0.14 - 0.07 fm explored in Ref. 0. 
We employ lattices with a physical extent of La ~ 3 fm with which we expect negligible 
finite-size effects: no significant effect is observed for La >2 fm within statistical error of 
about 2% [Q. We select five values of the hopping parameter k at which quark mass takes 
values corresponding to 171^,/ nip ^ 0.75, 0.7, 0.6, 0.5 and 0.4, the last point being closer 
to the chiral limit than previous attempts with the Wilson action. We abbreviate the two 
heavy quarks and the three light quarks as Si and Ui, respectively. 

Gauge configurations are generated by the 5-hit pseudo-heat bath algorithm with sub- 
sequent four over-relaxation updates. Quark propagators are solved on configurations fixed 
to the Coulomb gauge, using both the point source and an exponentially smeared source 
of a radius vq ^ 0.3 fm. For extraction of hadron masses, we use hadron propagators con- 
structed from only smeared sources, since they exhibit the earliest onset of effective mass 
plateaus and the smallest statistical errors among various source combinations. Masses are 
determined by correlated fits for degenerate hadrons as well as non-degenerate ones of 
the type SiUj for mesons and SiSiUj and SiUjUj for baryons. 

Chiral extrapolation is the most delicate issue in controlling the systematic errors. For 
this extrapolation, most studies have employed low-order polynomials in quark masses. 
Quenched chiral perturbation theory (QxPT) |l7HTl|, on the other hand, predicts character- 
istic singularities in the chiral limit. 

For pseudoscalar (PS) mesons made of quarks of masses mi and 1712, QxPT mass formula 
reads 

^\s.i2= + "^2){1 - 5[ln(2y4mi/A^) 

+7712/ {m2 - nil) ln(m2/mi)]} 

+S(mi + m2)^ + 0(m3), (1) 
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where the term proportional to 5 represents the quenching effect. The following two quan- 
tities 

_ 2mi mls^^2 2m2 ^Pg,i2 



x = 2 + !^^l±^ln(I^), (3) 



are then related by?/ = l + 5-a; + 0{vp?\ 

Our PS meson mass data converted to x and y are shown in Fig. |l[ The quark mass 
determined from the axial Ward identity [12,131 {j^^^^ is used since it has neither quenched 



chiral logarithms |jT^ nor ambiguity due to chiral extrapolations. We observe a clear non- 
zero slope in the range b ~ 0.08 - 0.12. The ratio ?^P5 12 /"^ps' n used in Ref. P] receives a 
correction from the O(m^) term in (|I]), and the plot does not fall on a common line. 

Similarly, an examination of the relation y = 1 — (5/2) ■ x expected for the decay constant 
ratio y = fi2/{fnf22) gives a non-zero value of 5 ~ 0.08 - 0.16. 

Finally, we make a correlated fit of the PS mass results to the QxPT mass formula. For 
this fit the vector Ward identity quark mass m^^^ = (1/k — 1/kc)/2 is used, taking Kc as one 
of the fit parameters. The fitting, made independently for degenerate and non-degenerate 
cases, yields mutually consistent results for 6 in the range 6 ~ 0.06 - 0.12 when is varied 
over ^ 0.6 - 1.4 GeV. 

A noticeable result from this fit is that the value of Kc agrees well with n-^^^ at which lin- 
early extrapolated m^^^ vanishes with at most 2. 80" discrepancy, whereas quadratic (cubic) 
extrapolation of in m^^^ results in Kc differing as much as 17a (12(j). This observation 
suggests that the QxPT fit is legitimate and rriq^^ should vanish at the same k 

value from their definitions. 

We take these results as strong support that our PS meson results are consistent with 
the presence of quenched chiral logarithms with a magnitude 5 ~ 0.1. We should mention 
that there is in (|I|) another quenched singularity term of the form a$m^ Inrrij Our mass 
data do not show clear indications for the presence of this term. 

QxPT mass formulae for vector mesons 0] and baryons [|ll[] can be schematically written 



as 

mH{mps) = mo + Ci/2mps + Ciuips + C^/2m\g, (4) 

where Cj are polynomials of the couplings of the quenched chiral Lagrangian. The coefficient 
Ci/2, which is proportional to 5, represents the quenched singular term. Encouraged by the 
results for the PS mesons, we attempt to fit our mass data to @). Fully constraining the 
large number of coupling parameters (6 for vector mesons and 11 for baryons in addition to 
5 and «$) is difficult under the limitation of the accuracy of our mass data. We therefore 
employ 5 = 0.1 and = 0, and drop the couplings of the fiavor-singlet PS meson to vector 
mesons and baryons. Fits are made to degenerate and non-degenerate data together ignoring 
correlations, as we find correlated fits to be quite unstable. Errors are calculated with the 
single elimination jackknife method. 

With this procedure we obtain a good fit with small x^/dof < 0.8 keeping all 0{mps) 
and 0{m'j}g) terms for vector mesons and decuplet baryons. For octet baryons, we also 
include 0{m?pg) terms since the nucleon mass shows a clear negative curvature, which is 
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opposite to the predicted negative sign of C1/2 (see Fig. |^). We drop the octet-decuplet 
couphng because the decuplet baryon fits lead to a value of this coupling consistent with 
zero. These fits yield Cxi2 = -0.071(8) for p, -0.118(4) for nucleon and -0.14(1) for A 
when averaged over the four /3 values of the simulation. These values are much smaller than 
phenomenological estimates: C1/2 = —0.71 for p [|1^] and —0.27 for nucleon ||11|| . 

We have also attempted to include the 0{m?pg) terms in the fit for vector mesons and 
decuplet baryons, and the octet-decuplet coupling for octet baryons, since it is theoretically 
more consistent to include them. Making this more general fit, we find the values of couplings 
to change significantly. In particular the coefficient of the 0{mps) term becomes C1/2 ~ 
—0.29 for p, —0.01 for nucleon and —0.02 for A. The stability of fit worsens considerably, 
however, compared to the fit above: the covariance matrix of the fit is close to singular, and 
some of the couplings exhibit significant variation with the lattice spacing. Nonetheless, as 
we illustrate in Fig. |^ for representative hadron masses, the curves from the two fits (dashed 
lines for the general fit and solid lines for the simplified fit as described in the previous 
paragraph) are indistinguishable in the range of measured points. The deviation remains 
small even at the physical point (second points from left in Fig. |^). 

We conclude that the accuracy of the mass results and the covered range of quark mass 
are not sufficient to establish the presence of the quenched singular term for vector mesons 
and baryons. We adopt the former fit with simplified QxPT formulae in the subsequent 
analyses as it is overall more stable. 

We emphasize that adopting the latter general fit does not affect our conclusions below, as 
changes in hadron masses at the physical point are 5% (5cr) or less at finite lattice spacings 
and at most 1.2% (l.Sa) after the continuum extrapolation. We have also checked that 
varying 6 over 6 ~ 0.08 - 0.12 («$ over ^ —0.7 — 1-0.7) changes the mass results by no 
more than 0.4% (1.3a) (2.9%(4.7or)) at finite lattice spacings and 0.3% (0.3o-) (2.2%(1.4(t)) 
in the continuum limit. 

The physical point for up and down quarks and the lattice scale are fixed from the 
experimental values of = 0.1350 GeV and nip = 0.7684 GeV, and the strange quark 
mass by that of mx = 0.4977 GeV or = 1.0194 GeV. We then extrapolate the results 
linearly in the lattice spacing to the continuum limit. We find that x^/dof of the fit with 
only the leading scaling violation term is reasonably small (< 1.6) for each hadron. The 
statistical error in the continuum limit is 1-3%. In Fig. ^ are shown typical continuum 
extrapolations. The magnitude of scaling violation at a = 0.075 fm, the mid-value of the 
range of our lattice spacing, is at most 10%. We then expect that O(a^) terms have an effect 
of 1% or less, which is smaller than the statistical error. 

In order to examine how much our results in the continuum limit depend on the QxPT 
fit, we repeat the analysis with the conventional polynomial fit of hadron masses in terms 
of 1/k. This fit (cubic for nucleon, quadratic for others including PS mesons) yields masses 
at the physical point differing by at most 3% at finite lattice spacings relative to the QxPT 
fit, and less than 1.6% (1.6(t) in the continuum limit. The continuum extrapolations for the 
QxPT and polynomial chiral fits are compared in Fig. ^ with filled circles and solid lines 
for the former and open circles and dashed lines for the latter. 

We present our final result for the quenched light hadron spectrum in Fig. ^. The 
numerical values are given in Table Systematic errors, being comparable or less than the 
statistical ones, are not included. 
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Our chief finding, quite clear from Fig. H, is that the quenched spectrum systematically 
deviates from the experimental spectrum. Quantitatively stated, if one uses the K meson 
mass to fix the strange quark mass, the vector meson masses ttik* and are smaller by 4% 
(4(t) and 6% (5ct), respectively, the octet baryon masses are smaller by 6-9% (4-7(t), and so 
are the decuplet mass splittings (30% on average). Alternatively, if is employed to fix 
the strange quark mass, niK* agrees with experiment within 0.8% {2a) and the discrepancies 
for baryon masses are much reduced. However, rriK is larger by 11% (6cr). As one sees in 
Table |n| this 11% represents the largest deviation between our results and the experiment. 

Our finding of a small mass splitting between K and K* differs from that of Butler et al. 
who found an agreement with experiment. Another difference is the masses of H* and Vl 
baryon determined with irtK as input, for which they reported values higher than experiment 
by 3-5%, while our results are smaller by a similar magnitude. How these differences arise 
are shown in Fig. ^ where the results of Ref. are plotted by open triangles. For the 
nucleon both the results from Butler et al. [Q] and the MILC collaboration are consistent 
with experiment; our value is smaller by 7% (2. Sex). 

We also calculate PS decay constants and quark masses using tadpole-improved one-loop 
values for renormalization constants. For the PS decay constant we find /tt = 120.0(5.7) MeV 
and fx = 138.8(4.4) MeV with m^- as input, which are smaller than experiment by 9%(2cr) 
and 13%(5(t), respectively. Quark masses are determined by a combined linear continuum 
extrapolation of rn^^^ and m^^^, since the large difference of values from the two defini- 
tions at finite lattice spacings |[T3|JT5[] vanishes toward the continuum limit ||^. We obtain 
m„,rf=4.55(18) MeV, and m,=115.1(2.3) MeV {ttik input) or 142.9(5.8) MeV (m^ input) in 
the MS scheme at = 2 GeV. A 20% disagreement between the two values for originates 
from the small meson hyperfine splitting, and hence represents a quenching effect. 

In conclusion we have found that the light hadron spectrum in quenched QCD system- 
atically deviates from the experimental spectrum when examined with accuracy better than 
a 10% level. In the course of analyses we have observed a strong support for the presence of 
quenched chiral singularities for pseudoscalar mesons. Whether vector mesons and baryons 
also have such singularities, however, remains as a problem for future investigations. 

We thank all the members of the CP-PACS Project with whom the CP-PACS computer 
has been developed. Valuable discussions with M. Golterman and S. Sharpe are gratefully 
acknowledged. This work is supported in part by the Grants-in-Aid of Ministry of Education 
(Nos. 08NP0101 and 09304029). GB, SE, and KN are JSPS Research Fellows. HPS is 
supported by JSPS Research for Future Program. 
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FIGURES 
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FIG. 1. Results of a test for the presence of the quenched chiral logarithms in pseudoscalar 
meson mass data. See text for details. 
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FIG. 2. Degenerate hadron masses vs. m'pg at /3=5.9. The leftmost points arc values extrap- 
olated to the chiral limit, and the second ones from left are those at the physical point. Fitting 
curves from two types of chiral extrapolations based on QxPT are reproduced. See text for details. 
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FIG. 3. Typical continuum extrapolations of masses with mx as input. Filled (open) circles 
are data from the QxPT (polynomial) fits. Open triangles are from Butler et al. Q, the leftmost 
ones being their final results after finite size corrections. Stars represent experimental values. 




FIG. 4. Final spectrum results compared to experiment. 
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TABLES 



TABLE 1. Simulation parameters. 



P L^T a-^(GeV) La(fm) #conf. sweep/conf 

K90 3^56 1.934(16) 3.26(3) 800 200 

6.10 40^ 70 2.540(22) 3.10(3) 600 400 

6.25 48^ 84 3.071(34) 3.08(3) 420 1000 

6.47 64^112 3.961(79) 3.18(6) 150 2000 



TABLE II. Spectrum results. Deviation from experiment with its statistical significance is 
also given. 







rriK input 


input 






exp. 


mass(GeV) 


deviation 


mass(GeV) 


deviation 


K 


0.4977 






0.553(10) 


11.2% 5.6(7 


K* 


0.8961 


0.858(09) 


-4.2% 4.3a 


0.889(03) 


-0.8% 2.3(7 




1.0194 


0.957(13) 


-6.1% 4.8c7 








0.9396 


0.878(25) 


-6.6% 2.5(7 


0.878(25) 


-6.6% 2.5(7 


A 


1.1157 


1.019(20) 


-8.6% 4.7(7 


1.060(13) 


-5.0% 4.1(7 




1.1926 


1.117(19) 


-6.4% 4.1(7 


1.176(11) 


-1.4% 1.5(7 




1.3149 


1.201(17) 


-8.7% 6.8(7 


1.288(08) 


-2.0% 3.5(7 


A 


1.2320 


1.257(35) 


2.0% 0.7(7 


1.257(35) 


2.0% 0.7(7 


S* 


1.3837 


1.359(29) 


-1.8% 0.9(7 


1.388(24) 


0.3% 0.2(7 


"* 


1.5318 


1.459(26) 


-4.7% 2.8(7 


1.517(16) 


-1.0% 0.9(7 




1.6725 


1.561(2-1) 


-6.7%: 1.7(T 


1.617(10) 


-1.5% 2.6(7 
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